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INTRODUCTION 


This  presentation  offers  a mathematical  framework  for  a series  of  problems 
normally  encountered  under  the  heading  "sound  ranging,"  and  developments 
are  directed  toward  this  end.  However,  the  techniques  may  be  subject  to 
a broader  range  of  interpretation.  The  language  is  intended  to  be  strongly 
suggestive  of  actual  application  and  will  be  kept  of  such  a nature  as  to 
agree  with  intuitive  "feel."  As  considered  here,  the  fundamental  problem 
of  "sound  ranging"  will  be  the  location  of  a "sound  source"  when  the  ap- 
proximate position  of  an  expanding  wave  front  is  known  at  a number  of  times; 
or  more  appropriately , if  the  approximate  time  of  arrival  of  such  a wave 
front  is  known  at  a series  of  prescribed  locations  in  space. 


The  approach  taken  is  nonlinear  and  consists  of  two  distinct  steps:  a 
preliminary  estimation  of  sound  source  and  influencing  parameters,  and 
an  iterative  technique  for  refinement  of  such  estimates. 


Of  particular  Tntere&t_Js  the  interplay  which  arises  between  points  cho- 

sen  as  "receivih<T*tations,"  Jound  source  position,  accuracy  of  arrival 

’*'•••  1 f 

times,  and  integrity  ^n^^jMrologi cal  information. 


The  guidelines  choserTfor  tfie  mathematical  framework  were  sufficiently 

, j.  I 

general  to"a)j|jp»y(for consideration  of  many  differing  sound  wave 

v,  7 

fronts  without  destroying  an  undue  amount  of  intuitive  flavor.  Several 
1nt^h?*ting  and  occasionally  surprising  relationships  emerge  which  are 
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not  generally  evident,  especially  with  regard  to  the  role  of  microphone 
placement. 

The  main  flow  proceeds  from  the  qeneral  to  the  specific,  with  an  ultimate 
goal  the  improvement  of  overall  "sound  ranging  capability."  This  goal 
should  be  approached  by  the  augmentation  of  current  or  proposed  techniques 
with  a class  of  "vernier"  methods.  The  authors  believe  that  these  techniques 
may  be  initiated  at  a modest  computational  level  and  will  be  consistent 
with  modern  concepts  of  automation. 

In  this  initial  presentation,  a general  mathematical  skeleton  is  devel- 
oped and  then  applied  in  a simple  problem.  The  example  presented  will 
consider  a spherical  sound  wave  propagating  in  a medium  of  uniform 
temperature  and  in  the  presence  of  a wind  of  fixed  magnitude  and  direc- 
tion. This  model  will  be  sufficient  to  demonstrate  several  pertinent 
features,  yet  the  lack  of  complication  provides  ease  of  presentation 
and  interpretation.  As  more  complex  models  are  examined,  they  may  be 
contrasted  in  light  of  their  differences  from  the  basic  model  and  per- 
haps in  this  manner  be  more  easily  interpreted  themselves. 

To  this  point,  random  variations  within  the  structure  have  not 
been  considered.  However,  the  authors  have  explored  this  area  to  the 
extent  that  a probabilistic  approach  is  expected  to  be  highly  rewarding, 
both  from  the  point  of  view  of  sensitivity  analysis  and  from  the  stand- 
point of  information  recovery  versus  receiver  accuracy.  The  authors 
believe  stochastic  modeling  will  provide  a fruitful  avenue  of  approach, 
and  preliminary  results  are  encouraging. 
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THE  SOUND  RANGING  PROBLEM,  A GENERAL  MATHEMATICAL  APPROACH 
Fundamental  in  the  structure  to  be  developed  is  the  notion  of  an  extended 
wave  generating  function.  For  the  highly  specialized  application  encoun- 
tered here,  the  following  definitions  will  be  employed. 

Definition  1 

Let  w be  a real -valued  function  whose  domain  is  a compact  connected  re- 
gion of  Euclidean  4-space  ( E4 ) and  whose  range  contains  the  real  number 
zero.  Assume,  moreover,  that  w is  of  class  C2,  that  is,  possesses  con- 
tinuous partial  derivatives  of  the  second  order.  Then  w will  be  called 
an  extended  wave  generating  function. 

Let  V denote  the  domain  of  definition  of  an  extended  wave  generating 
function  w.  A typical  point  in  V will  be  written  as  (x,y,z,t);  with  D = 
{(x,y,z)  | (x,y,z,t)  e V}  carrying  a connotation  of  "points  in  space," 
while  T = {t  | (x,y,z,t)  e V } has  a context  of  "time."  The  notation  for 
V and  the  projections  D and  T will  be  consistently  utilized  throughout 
the  discussions  to  follow. 

Observe  that  the  continuity  of  projection  maps  assures  that  both  T and  D 
are  compact  and  connected;  hence  D may  be  visualized  as  a bounded  volume 
in  Euclidean  3-space  (E3),  and  T as  a real  interval. 

Attention  will  now  be  focused  on  sets  of  the  form 

W = {(x,y,z)  | (x,y,z,t)  e V;  w(x,y,z,t)  = 0}.  Some  descriptive  termi- 
nology concerning  these  sets  is  given  in  the  following  definition. 
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Definition  2 


Let  w be  an  extended  wave  generating  function  defined  on  the  region  V of 
E4.  For  each  t e T,  the  set  of  points  in  E3  defined  by  Wt  = 

{(x.y.z)  | (x,y,z,t)  e V ; w(x.y.z.t)  = 0}  will  be  called  the  wave  front 
at  t generated  by  w.  A point  (x.y.z)  in  Wt  will  be  said  to  lie  on  the 
appropriate  wave  front. 


Suppose  that  {w^  | j e J}  is  an  indexed  family  of  extended  wave  generating 
functions;  and  let  {Dj  | j e J},  {D.  | j e J},  and  (T\  | j e J}  have  their 


obvious  interpretations.  Let  (x.y.z)  be  a point  in  D.,  which  will 

j c j j 

be  assumed  to  be  nonempty.  If  there  exists  (t.  | t.  e T.;  j e J}  such 
S~\  J •'  J 

that  (x.y.z)  c j c j(Wj)t  , where  for  each  .i  e J,  (Wj)t  denotes  the 

wave  front  at  tj  generated  by  w^,  the  point  (x.y.z)  will  be  called  a 

point  of  contingency  for  the  family  (w.  | j e J}.  A special  point  of 

J 

this  nature  will  be  of  interest  in  the  solution  of  the  sound  ranging 
problem  as  presented  in  the  current  context. 


At  this  point,  it  will  be  expedient  to  define  a special  point  which  will 
emerge  as  a point  of  contingency  for  certain  families  of  extended  wave 
generating  functions.  These  will  arise  in  a logical  manner  from  a single 
wave  generating  function,  and  a known  family  of  discrete  points  in  the 
domain  of  this  latter  function.  This  point  will  be  called  a source 
point,  and  its  description  is  the  content  of  the  following  definition. 


ir 
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Definition  3 


Let  w be  an  extended  wave  generating  function  defined  on  V.  For  any 
(x,y,z)  e D,  let  r(x,y,z)  designate  the  set  of  all  linear  graphs  in  E3 
with  origin  (x,y,z). 

Now  suppose  that  (xQ,  yQ,  zQ,  tQ)  is  a point  in  V with  the  following 
properties 

a.  For  any  y e r,  t e T,  y 

b.  For  any  y e r,  and  sequence  (tn  | tfl  e T;n  1,2,...}  such  that 

tn  = tQ;  if  (pn,qn,rn)  is  in  >PlWt  , it  follows  that 

'r  II  <Wrn>  • <VVzo>  11-0. 

Then  (x0»y0»z0»t0)  will  be  called  a source  point  of  w. 

Theorem  1.  The  source  point  as  designated  by  Definition  3 is  unique; 

that  is,  an  extended  wave  generating  function  can  have  at  most  one  source 

point. 

Proof.  Suppose  there  exist  two  distinct  source  points  (x1 ,y]  ,z]  ,t] ) , 
( x2 ,y2»z2*^*2^ ‘ For  times  sufficiently  near  t1  and  t2,  all  wave  fronts 
belonging  to  these  times  may  be  separated  by  open  spheres.  By  choosing 
a point  outside  either  of  these  spheres  so  that  a line  connecting  this 
point  to  the  center  of  one  of  the  spheres  is  not  tangent  to  the  other, 
part  a of  Definition  3 Is  violated.  This  proves  the  theorem  and  demon- 
strates uniqueness  of  the  source  point. 


A fundamental  property  of  extended  wave  generating  functions  with  source 
points  which  will  be  of  considerable  importance  in  the  development  to 
follow,  and  which  will  offer  one  approach  to  the  problem  of  numerical 
determination  of  source  points,  can  now  be  demonstrated. 

Theorem  2.  Let  w be  a wave  generating  function  with  source  point 
(Vyo>zo>t0).  Let  J be  a finite  index  set,  and  for  j e J,  let 
(x-,y ..z-.tj  lie  in  V.  Assume  that  {(x^.y-.z-.t.)  | j e J}  forms  a 

J J J J \J  \J  SJ  \J 

distinct  set  of  points.  It  is  now  possible  to  define  a family  of  trans- 
formations from  E4  onto  E4  by  letting  tj  for  each  j e J be 


TJ:<^ 


■x'  = X - (Xj  - xQ) 

y’  = y - (yj  - y0) 

2'  ■ 2 - (zj  - v 


Lt'  -‘-w 

For  each  j e J,  t.  possesses  a differential  equal  to  the  appropriate 

J 

identity  matrix;  therefore,  the  Jacobian  is  of  constant  value  1 through 
E4.  In  particular,  is  a pure  translation  in  E4,  and  hence  carries 
compact,  connected  sets  into  sets  of  similar  nature.  Now,  for  each 
j e J,  let  w.  be  defined  for  each  p'  e t.(P)  by  w.(p')  = w(t_1(p')). 

J J J 

Clearly  is  well-defined  for  each  j e J.  It  follows  that  {w^  | j r J} 
is  a family  of  extended  wave  generating  functions,  and  (x0>y0»z0)  is  a 
point  of  contingency  for  this  family. 
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4 Proof.  As  has  been  seen,  the  nature  of  the  family  of  translations 

{t.  | j c J}  assures  that  each  member  of  (w.  | j e J}  is  a real-valued 

J J 

function  defined  on  a compact  connected  region  of  E4.  It  also  follows 
that  each  is  of  class  C2,  since  w possesses  this  nature. 

Observe  now  that  for  any  j t J,  (x.,y.,z.,t .)  = p.  is  mapped  by  t.  onto 

J J J J J J 

(xo’yo’zo,to)  ' V hence’  P0  is  in  the  domain  of  Wj  for  each  j e In 
addition,  Wj(pQ)  = w(t_1(po))  = w(p^)  = 0;  therefore,  zero  is  in  the 

range  of  w^.  for  each  j e J and  each  w^  is  an  extended  wave  generating 
function.  Examination  also  yields  that  Cx0»y0»z0)  lies  in  (Wj)t  ; that 
is,  the  wave  front  at  t.  generated  by  w.;  and  this  for  each  j e J.  It 

J J 

therefore  follows  that  (.xQ,yQ,zQ)  is  a point  of  contingency  for 
{w.  | j e J}.  This  proves  the  theorem. 

i| 

In  reality,  an  important  corollary  has  essentially  been  proved. 

Corollary  2.1.  If  w,  pQ,  and  (w^  | j e J}  are  as  given  in  Theorem  1, 
the  source  point  of  w lies  among  the  common  zeroes  of  the  family 
(Wj  | j e J}. 

Corollary  2.1  links  the  problem  of  determining  a source  point  of  an  ex- 
tended wave  generating  function  to  the  problem  of  determining  the  proper 

zero  of  (Wj  | j e J}.  As  will  be  seen  later,  a series  of  computational 
devices  for  determining  source  points  may  arise  from  this  fact. 

At  this  point,  a precise  definition  of  the  sound  ranging  problem  may  be 
phrased  in  current  context. 
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Definition  4 


Let  w be  an  extended  wave  generating  function  and  {p.  | j e 0}  be  a 
finite  family  of  points  in  V with  {(Xj.yj.Zj)  | j e J)  distinct  in  E3. 

The  problem  of  determining  the  source  point  of  w,  if  such  exists,  will 
be  called  a sound  ranging  problem. 

In  the  study  of  general  sound  ranging  problems,  an  especially  revealing 
form  of  the  extended  wave  generating  function  is  described  in  the  following 
definition. 


Definition  5 

Let  w be  an  extended  wave  generating  function  having  the  following 
properties 

a.  For  any  t,t'  in  T,  Wt  f Wt, 

b.  There  exists  (x,y ,z)  e D having  the  property  that  if  y e r (x,y,z), 
and  t,  t'  are  in  T with  t <_  t ' ; and  if  in  addition,  (p,q,r)  e y 

Wt  and  (p'.q'.r1)  e y/^Av^,,  then 

||  (x,v,z)  - (p,q,r)  II  < | | (x,y,z)  - (p'.q’.r')  ||. 

Then  w is  said  to  be  an  extended  wave  generating  function  with  wave  fronts 
expanding  about  (x,y,z).  If  strict  inequalities  hold  in  b the  wave  fronts 
in  question  will  be  said  to  be  strictly  expanding. 

The  following  theorem  exhibits  some  of  the  more  salient  features  of  ex- 
tended wave  generating  functions  with  wave  fronts  expanding  about  a source 
point. 
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Theorem  3.  Let  w be  an  extended  wave  generating  function  with  source 
point  pQ,  and  which  generates  wave  fronts  expanding  about  (x0»y0»z0)- 
Then  the  following  conditions  hold 

a.  The  set  W^.  must  contain  only  the  point  (xQ,yo,zo). 

b.  The  point  tQ  is  the  greatest  lower  bound  of  T. 

In  this  case,  tQ  is  called  the  blast  time;  and  (x0«y0»z0)  the  point  of 
origin  of  the  expanding  family  of  wave  fronts. 

Proof:  Let  t'  be  given  such  that  t1  £ tQ.  As  w generates  wave 

fronts  expanding  about  (x0»y0»z0)»  it  must  follow  tnat  ||  (xQ,yo,zo)  - 
(P.q.r)  ||  = 0 for  any  (p,q,r)  e Ht,,  hence  Wt,  h {(xo,yo>zo)}.  By  con- 
dition a of  Definition  4,  it  must  follow  that  t'  = t . Conditions  a and 
b follow  at  once. 

It  might  be  recalled  from  Corollary  2.1  that  a source  point  of  an  extended 
wave  generating  function  lies  among  the  conrnon  zeroes  of  a family  of  func- 
tions of  class  C2,  each  defined  on  a compact,  connected  region  of  E4.  If 
the  entire  set  of  zeroes  of  this  family  may  be  determined,  there  still 

remains  the  problem  of  isolating  the  source  point  from  the  remaining 
zeroes.  If  a set  of  restraints  is  known  which  assures  this  separation, 
the  resulting  system  will  produce  what  will  be  called  a completely 
determinate  sound  ranging  problem.  As  will  be  seen,  however,  the  prob- 
lem of  selecting  a proper  zero  may  be  approached  from  another  avenue. 
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The  problem  of  solving  a set  of  simultaneous,  nonlinear  equations  of  the 
form  (w.j(x,y,z,t)  = 0 | j e J}  is  well  known  in  the  literature  of  numerical 
analysis,  and  an  introductory  presentation  is  given  for  example  in  [2]. 

The  particular  approach  which  will  be  examined  in  context  with  the  pre- 
viously stated  problems  will  be  the  familiar  Newton-Raphson  iterative 
technique,  or  a slight  variant  thereof,  utilizing  the  concept  of  the 
Moore-Penrose  Pseudoinverse  [3]. 

IP 

The  notation  to  be  followed  hereafter  is  largely  that  of  [1],  and  has  the 
advantage  of  smoothness  in  manipulation  and  interpretation.  In  this  nota- 
tion, if  (wi  | i = l,2,....n}  is  a family  of  real -valued  functions  defined 
on  a region  of  E4,  and  of  class  C1  or  greater,  this  defines  a transforma- 
tion, call  it  T,  from  E4  into  En.  For  each  p in  the  comnon  domain  of  the 
family  (w^  | i = 1,2,... .n},  an  n x 4 matrix,  cal  led  the  differential  of 
the  transformation  (in  symbols  dT(p)),  is  defined  by  letting  row  k be 

3wk  3wk  9wk  3wk 

given  by  ~ | p,  ~ - | p,  | p,  | p,  for  each  k among  1,2, n. 

In  other  words,  dT(p)  is  the  matrix  whose  kth  row  is  the  qradient  of  w^ 
at  p. 

If  T is  as  previously  defined,  but  with  n = 4,  and  p‘  is  any  zero  of  T; 
that  is,  T(p')  is  the  zero  vector  in  E4,  there  is  a neighborhood  of  p* 

such  that  for  any  pj  in  this  neighborhood,  the  sequence  (pk  | k = 1,2, } 

defined  for  k = 1,2,...,  by 
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"ktl  = "k  - [dT(Pk)]''  T(Pk> 


0) 


will  converge  to  p'  if  the  determinant  of  dT(p')  (in  symbols  | dT(p')  | ) 
is  nonzero. 

A more  general  scheme,  and  one  which  will  prove  of  wider  application  than 
the  above,  may  be  defined  by  utilizing  the  concept  of  the  Moore-Penrose 
Pseudoinverse.  This  allows  n to  be  any  finite  number  at  least  four. 
Advantages  of  such  a scheme  are  obvious,  but  caution  must  be  exercised  in 
the  numerical  inversion  of  the  involved  matrix;  as  a worsening  of  the 
condition  of  the  linear  system  may  occur  as  n becomes  large.  Such  a 
situation  may  be  accompanied  by  computational  difficulties.  This  scheme 
will  be  discussed  more  fully  as  matters  progress. 

WAVE  GENERATING  FUNCTIONS  WITH  ONE  OR  MORE  PARAMETERS 
Let  n(V)  be  the  set  of  all  wave  generating  functions  on  E4  having  a com- 
mon domain  V.  Let  A be  a region  of  Euclidean  m-space  which  is  compact 

and  connected.  Interest  will  be  with  a function  f whose  domain  is  A and 
whose  range  is  fi(P).  The  set  f{A>  = {f(a)  | a e A};  that  is,  the  ranqe 
of  f will  be  called  the  family  of  wave  generating  functions  with  param- 
eters in  A.  If  a Is  in  A and  p c V,  the  value  of  f (a)  at  p will  be  de- 
noted by  [f(a)](p). 

The  set  f[A]  will  be  said  to  define  a family  of  wave  generating  functions 
with  continuous  parameters  if  for  any  a,  b in  A and  any  real  e > 0,  there 
exists  6 > 0 such  that  | [f(b)]  (p)  - [f(a)]  (p)  | < e when  ||  b - a ||  < 
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and  this  condition  holds  for  any  p e D.  In  a similar  fashion,  one  may  de- 
fine wave  generating  functions  with  absolutely  continuous  parameters,  of 
class  C1,  C2,  etc.,  the  interpretation  being  straightforward. 

A source  point  of  f (a)  may  not  be  a source  point  of  f(b)  if  a f h.  More 
attention  will  be  given  this  important  concept  later  in  this  discussion 
where  the  applications  to  the  sound  ranging  problem  will  be  clarified. 


AN  APPLICATION  OF  THE  MOORE-PENROSE  PSEUDOINVERSE 
Let  A be  a region  of  Em  and  let  f:A  n(P)  define  a family  of  wave  gen- 
erating functions  with  parameters  a = (a  ,a  ,. . .a  ) t A.  It  will  be 

l 2 m 

assumed  that  the  family  f{A}  is  of  class  C1  or  greater.  As  has  been 
discussed  previously,  any  member  of  f { A)  is  a wave  generating  function, 
and  for  a,  b In  A with  a / b,  the  source  point  of  f (a)  will  not  in 
general  be  equal  to  that  of  f(b),  and  it  is  possible  that  such  a point 

fails  to  exist  for  either.  However,  the  assumption  will  now  be  made  for 
each  at  A,  f (a)  has  a source  point. 


Now  suppose  for  some  at  A,  a family  of  n points  { (x-  ,y^,z^,t^)  | 

1 = 1,2, n)  exists  with  n >_  4.  Suppose  it  is  further  known  that 

Wj  Is  the  wave  front  at  t^  of  f (a),  (x^ ,y^,z^)  lies  in  , and  that 
,y.j  ,z.j  ,t.j ) Is  distinct.  One  may  then  define  a transformation  of 
Into  En  exactly  as  done  previously,  obtaining: 

‘ f^aHx'  ,y',z',t') 

f (aHx'.y'.z'.f) 

2 


T * < 


if 

each 

E- 


f n (a ) ( x ' ,y',z'  ,t') 
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Note  that  the  set  of  points  p such  that  T(p)  = 0 contains  the  source 
point  of  f(a)  which  will  be  called  pQ(a). 

It  is  not  overly  difficult  to  show  that  if  p is  chosen  within  a suf- 

l 

ficiently  small  neighborhood  of  pQ(a),  the  iterative  scheme  defined  by 

Pk+i  = Pk  - t[dT(pk)]C  [dT(pk)]}_1  [dT(pk)J*  T(pk)  (2) 

will  converge  to  pQ(a)  if  dr(pk)  is  always  of  rank  at  least  4.  This 
fact  will  be  of  interest  in  several  practical  applications  and  will  be 
of  particular  interest  when  a random  model  is  to  be  considered. 


At  this  time,  suppose  that  n >_ m + 4,  and  that  f (A)  is  of  class  C1  with 

respect  to  the  parameters  (a  ,a  a„)  e A.  Consider  the  transforma- 

12  n 

tion  from  (m  + 4)-space  to  n-space  defined  by 


t:  < 


'fi(a)(x,y,z,t) 

f (a)(x,y,z,t) 

2 


(3) 


[fnU)(x,y,z,t) 

The  Iterative  procedure  defined  by  choosing  p^  = (xj  ,zj *aj » 
a '.....a^)  and  for  k > 1 ,pk  by 


Pk+1  * Pk  - ([dT'(pk)]f  [dT’(pk)]}_1  [dT'(p^)]t  T'(pp  (4) 
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will  converge  to  = (x0*y0,zo,to,ai,a2*‘*‘  'am^  if  ^ pk  " po  ^ m+4 
Is  sufficiently  small  and  dT1  is  of  rank  m + 4 at  the  source  point  of 

f (a) . Note  what  this  states.  If  f{A}  is  a family  of  extended  wave 
generating  functions  with  parameters  (a^ ,a^ ,. .. .an)  e A,  one  may  deter- 
mine the  source  point  of  f(a)  and  the  parameters  (a^ ,a2 aR)  from 

knowledge  of  points  and  the  wave  fronts  of  adjoining  or  nearby  wave 
generating  functions.  In  a physical  situation,  it  means  that  if  suf- 
ficient devices  are  placed  in  the  field  to  determine  necessary  arrival 
times,  these  times  may  contain  information  about  parameter  sets  describ- 
ing the  family  of  extended  wave  generating  functions,  as  well  as  informa- 
tion concerning  the  source  point. 

Applications  of  the  Moore-Penrose  Pseudoinverse  will  prove  to  be  of 
special  interest  when  one  considers  various  random  models  for  field 
application,  and  vulnerability  is  a factor  to  be  considered.  It  might 
be  noted  that  in  the  case  that  n = 4,  expressions  (1)  and  (2)  coincide, 
and  if  n = (m  + 4),  expressions  (1)  and  (4)  coincide.  More  specific 
information  can  be  derived  if  a particular  example  is  considered  for 
f(A},  which  leads  into  the  following  section. 

SOUND  RANGING:  A SIMPLE  MODEL  AND  SOME  RAMIFICATIONS 
Consider  the  following  highly  idealized  situation.  A sound  source  is 
located  at  some  point  in  a region  of  space  over  which  is  specified  a 
rectangular  Cartesian  coordinate  system  (X,Y,Z).  It  will  be  assumed 
that  a steady  wind  of  velocity  (u,v,q)  is  blowing  throughout  the  region 
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of  interest,  and  that  the  temperature  is  a constant  k degrees  Kelvin. 

A candidate  function  for  describing  the  evolution  of  sound  wave  fronts 
has  been  suggested  by  some  authors,  for  example  [4],  as 

C(x  - X0)  - u(t  - tQ)]2  + [(y  - y0)  - v(t  - t0)p 

+ [(2  - zQ)  - q(t  - tQ)]2  - c2k(t  - tQ)2  = w(x,y,z,t)  (5) 

where  c is  of  constant  value. 

If  the  region  concerned  is  a compact  (closed  and  bounded)  connected 
region  of  E4;  w,  as  defined  by  (5),  is  an  extended  wave  generating  func- 
tion. In  fact  the  following  theorem  may  be  proved. 

Theorem  3.  Let  w be  as  given  by  condition  (5)  and  be  defined  on 
a compact  connected  region  of  E4  containing  the  point  (x0»y0»z0»t0)* 

Then  w is  an  extended  wave  generating  function  with  source  point 
(x0,y0,z0,t0).  If  /u2  + v*  + q2  <.  c that  is,  if  the  wtndspeed  is 
less  than  Mach  1;  and  in  addition,  tQ  is  the  greatest  lower  bound  of  T, 
then  w generates  a family  of  wave  fronts  expanding  about  (x0*y0»z0)'.  I f 
windspeed  is  less  than  Mach  1,  the  family  of  wave  fronts  is  strictly 
expanding. 

Proof.  Since  by  hypothesis,  w is  defined  on  a compact  connected 
4 

region  of  E , the  fact  that  w is  an  extended  wave  generating  function 
is  a matter  of  routine.  It  is  clear  from  the  geometry  that  any  linear 
trace  whose  origin  is  at  (x0*y0»z0)  must  intersect  the  spherical  surface 
wt  for  any  t c T such  that  t > tQ,  and  that  any  sequence  of  times  in  T 
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have  property  b of  Definition  3.  Hence  (x0*y0»z0*t0)  is  the  source 
point  of  w.  The  remaining  verifications  are  equally  straightforward  and 
are  left  to  the  reader. 


This  extended  wave  generating  function  may  be  visualized  as  giving  rise 
to  a family  of  undeformed,  concentric  spherical  wave  fronts  propagating 
outward  at  a rate  of  c A,  and  being  translated  along  (u,v,q)  with  a 
speed  of  /u2  + v2  + q2. 


As  has  been  mentioned,  this  is  not  a complex,  or  even  a particularly 
realistic  "real  world,"  model;  and  many  physical  phenomena  which  are  nor- 
mally encountered  are  not  included.  However,  there  remains  an  adequate 
amount  of  Information  to  warrant  Its  careful  examination.  In  addition, 
in  such  an  uncomplicated  model,  one  is  not  overly  hampered  by  problems 

of  notation;  and  some  degree  of  intuitive  flavor  is  retained  which  could 
become  lost  if  a more  complicated  model  were  chosen. 


For  the  purposes  of  this  presentation,  all  quantities  will  be  assumed 
deterministic,  that  is,  not  random  in  nature.  The  ramifications  of  prob- 
abilistic models  of  varying  complexity  are  of  considerable  interest,  and 
sufficient  progress  has  already  been  made  to  assure  that  this  is  a highly 
productive  area,  with  certain  elegant  relationships  emerging.  However, 
sufficient  results  are  available  from  the  deterministic  model  to  bring  a 
number  of  significant  facts  to  the  surface. 


With  w as  given  by  relation  (5),  It  can  be  seen  that  for 
{ (xj ,yj ,Zj ,tj)  | j = 1,2, n},  a set  of  points  in  V with 
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{(x.,yi,zi)  | j = 1,2, n>  distinct,  one  may  define  a family 

J J J 

(Wj  | j = 1,2, n)  as  in  Theorem  2 by  utilizing  the  transformations 


rx'  = x - (Xj  - xQ) 


1 $2  9 • • • • n 


y'  = y - (yj  - y0) 

z'  = z - - z0) 

Lt’  ■ 1 - "j  - V 

Substitution  in  expression  (5)  yields,  for  j = 1,2, n 

Wj(x'  ,y*  ,z\t')  = wC-rjKx'.y'.z'.t,)] 

= Ex'  + x j - 2xQ  - u(t ' + t.  - 2tQ)]2 

+ [y*  + yj  - 2y0  - u(t'  + tj  - 2t0)]2 

+ [z*  + z.  - 2zq  - f(f  + t.  - 2t0)]2 

- c2k[t'  + tj  - 2tQ]2 

Observe  that  for  any  j among  l,2,....n. 


"j'WVV  ■ [(xj  ■ xo>  • u(tj  - ‘o'l2 

* t(»,  - V ■ '(‘j  - VJ2 

* CUj  - *„>  - - 10)]2 

- c!k[tj  - t0]2 

= 0 

This  may  be  derived  through  algebraic  manipulation  bearing  in  mind  the 
fact  that  (xj»yj*zj)  lies  on  the  wave  front  of  w for  each  j among 
1 ,2,. . . .n. 
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Now  suppose  that  n > 4,  and  that  the  parameters  u,  v,  q,  k are  known. 
Consider  the  transformation  T from  E4  into  En  defined  by 

-w^x.y.z.t) 

w (x,y,z,t) 

2 

T:<  (6) 

.wn(x,y,z,t) 

The  differential  of  T at  the  point  pQ  = (x0»yo,Zo’V  is  an  n x 4 matrix 

whose  jth  row  is  given  by 

<«<p#»lfJ  ■ zc(xj  - x0)  - uttj  - t0n 

(dT('’o»2,J  = 2[(yj  ' V ' v(tj  ‘ ‘o>3 

Mr* »3,J  ■ 2t(2j  - zo>  - <l(tJ  ' ‘o’3 

(dT(po))4,j  ' * v(dT(po))2,J(tJ  ‘ V 

+ q(dT(p0))3jj  - c2k(tj  - t0)3 

Recall  that  the  sequence  {pn  | n = 1,2, } as  defined  by  (2)  will  con- 

verge to  pQ  if  the  matrix  dT (pQ)  is  of  rank  4,  and  if  p^  is  chosen 
within  a sufficiently  small  neighborhood  of  p . Consider  the  following 
chain  of  reasoning.  The  matrix  dT(pQ)  is  known  to  be  of  unaltered  rank 
If  (dT(pQ))4  j is  replaced  by 

(<T(P0»4,J  * (ir(po))4,j  * "(dT(Po>>l,J  * ,(dT(*’o))2,j 
+ q(dT(p0))3j  * -2c2k(tj  - tQ) 

for  each  j among  1,2 n.  Through  a similar  reasoning  process,  one  makes 

the  following  substitutions 
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WfMlJ  ' (dTlpo»l,j  ' ^dT(Po»4,j 

<dT<Po))2,j  ' (dTM2.j  - ?k<dT<'’o))«,J 

(dT(l’o»J,j  " (dT(Po»3.J  - ^k(dT<’’o))4,J 

Note  that  the  matrices  dT(pQ)  and  dT'(pQ)  have  identical  rank. 

Finally,  the  observation  can  be  made  that  if  dT"(pQ)  is  defined  as 
dT"(p0)  = -2nc2k(dT' (pQ))»  the  rank  of  dT(pQ) ) and  dT"(po)  are  also 

identical,  and  dT"(p0)  may  be  written  as 

<dT"0’o»l,j  ■ txj  - x„) 

1 (dT"(P0)>2,J  ■ <*j  - V 

j = 1 ,2, n 

<dT"<',o»3,j  = (zj  - zo> 

(dT"<Po»4,J  ’ “j  - to) 

Now  suppose  that  the  points  {(Xj.yj.Zj.tj)  | j = 1,2, n}  are  of  such 

a nature  that  constants  a,b,c,d  may  be  found  which  are  not  all  zero,  and 
for  which 

a(xj  - xQ)  + b(yj  - yQ)  + c(Zj  - zQ)  + d(tj  - tQ)  = 0.  (7) 

It  follows  that  the  matrix  dT(pQ)  must  be  of  column  rank  3 or  less,  hence 
no  iterative  scheme  similar  to  the  one  given  by  (2)  can  converge  to  p . 

A special  case  of  the  above  which  appeals  to  geometric  intuition  occurs 
when  d of  condition  (7)  is  known  to  be  zero.  For  this  case,  condition  (7) 

Is  equivalent  to  the  statement  that  all  points  {(Xj.y^.Zj)  | j = 1,2 n} 

lie  on  a single  straight  line  through  the  point  (x0,yQ,z0).  This  condi- 
tion will  be  quickly  recognized  by  those  with  sound  ranging  experience. 
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Considering  alternately  the  conditions  a = 0,  b = 0,  c = 0,  together  with 
d = 0,  it  can  be  shown  that  dT"(p0)  and  therefore  dT(pQ)  will  be  of  rank 
less  than  4,  if  a straight  line  connecting  the  proper  points  can  be 
constructed  in  the  appropriate  projection.  For  example,  if  a line  con- 
necting {(x.,y.)  | j = 1,2, n}  and  (x  ,y  1 can  be  constructed  in  the 

j j oo 

(X,Y)  plane,  and  so  forth. 


It  would  appear  that  among  all  of  the  various  possible  configurations  for 
{(xj,yj,Zj)  | j = 1,2, n},  that  one  which  permits  the  above  set  of 

points,  or  projections  thereof,  to  lie  on  a straight  line,  is  intrinsi- 
cally of  more  than  usual  risk  from  a standpoint  of  source  point  determina- 
tions and  Information  recovery  through  a scheme  such  as  given  by  iterative 
scheme  2. 


RECOVERY  OF  PARAMETERS  ASSOCIATED  WITH 
A GIVEN  WAVE  GENERATING  FUNCTION 

Consider  now  the  following  situation.  Let  A be  a given  compact  connected 
region  of  E4.  Let  a = (a  ,a  ,a  ,a  ) be  a point  in  A,  and  define  f (a): 

12  3 4 

A -*■  E1  by  choosing  p * (x,y,z,t)  e V,  and  let 


Cf(a)](p)  * [(x  - xQ)  - at(t  - tQ)]2 

♦ t(y  - y0)  - «2(t  - t0n2 

+ C(Z  - z0)  - a3(t  - t0)]2 


(8) 


- c2a  [t  - tj2 

4 0 
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Let  V be  the  domain  of  the  wave  generating  function  given  by  (5).  Then 
for  each  a c A,  f(a)  e (P) , so  that  f is  well  defined  on  the  set  of 
wave  generating  functions  with  domain  P.  The  family  f(A},  then,  is  a 
set  of  wave  generating  functions  in  P with  parameters  in  A.  In  fact, 
(5)  is  a special  case  of  (8)  with  A = (u,v,q,k).  It  is  seen  from  (8) 
that  f(A>  is  also  a family  of  wave  generating  functions  with  parameters 
of  class  C2  over  A. 

Now  let  a'  c A be  given,  and  {p.  | j = 1,2, n)  be  a set  of  points  in 

P.  One  may  define  a transformation 

rtfu^i^p') 


UfU'jyp') 

by  allowing  [f («x 1 ) 3 j = [f(a‘  JlCxj^p’ )]  for  j = 1,2 n. 

Suppose  that  n >_  8.  If  (a^.p  ) is  in  a sufficiently  small  neighborhood 
°f  (a0,p0),  the  iterative  scheme  specified  by  (4)  will  converge  to 
(a0,P0)  if  the  rank  of  d H’(p0)  is  8 or  greater.  To  determine  con- 
ditions which  may  present  problems  with  regards  to  convergence  of  the 
iterative  scheme  mentioned  earlier,  procedure  will  be  as  follows,  v 
may  clearly  be  considered  a transformation  of  E8  into  En  of  class  C2. 
The  differential  of  this  transformation  at  UQ.P0)  may  be  calculated 
from  (8)  as  in  the  previous  section,  obtaining  columnwise 


22 


[d  Y(a'  ,pQ) 3-j  j = 2[(Xj  - xQ)  - - tQ)] 

[d  *(<t,.p0>]2iJ  ■ 2[(yj  - y0)  - - t0)] 

[d  P0)]3jj  = 2t(Zj  - a0)  - - t0)] 

[d  y(a,,po)]4>j  » -{2  a][d  «Ka'  ,pQ)]1  . + a£[d  *(a'  .P0)32  j 

* a3[d  * c2a4[tj  • lo],) 

[d  ’U',P0!]5iJ  • - (tj  - t0)[d  *(a'.P0)]M 

[d  *(<■•, p0n6ij  = - (^  - yu  *(«'.p0)]2>j 

[d  *U',P0)]7jJ  = - (tj  - t0)[d  t(d',p0)]3J 

[d  »(«•  .p0n8ij  * - t2(tj  - 10) 

for  j = 1,2 n. 

Observe  that  by  reasons  similar  to  those  employed  in  the  pre- 
ceding section,  the  matrix  d f(a',p0)  will  have  rank  the  same  as  a 
matrix  in  which  the  substitution 

[d  *(«' .p0)]4j  = [d  *(u'  ,p0>34  j * aj[d  »(a',p0)]1_j 

♦ aj[d  »(«',p0)]2fj  a aj[d  »<«• ,p0)]3j 
■ - c2a4<‘j  - V 

has  occurred. 
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One  may,  after  a series  of  substitutions  like  those  of 

the  preceding  section,  arrive  at  a matrix  d H,"(a',p0)  having  a rank 

Identical  to  d f(a' ,pQ)  but  with  components  given  columnwise  by 

td  »“(«'. »>„>],  ,j  * (Xj  - x0) 

[d  »"(»'. P0)]2J  * (yj  - y0) 

[d  *"(<■'. p0)]3iJ  * (Xj  - x0) 

[d  *■(«', P0)]4,j  * (tj  - t0) 

Cd  »”(»'. P0)]5,j  • (tj  - t0)(Xj  - x0) 

[d  P#)]6jj  ■ (tj  - t0)(yj  - y0) 

[d  y"(a'.Po)]7>J  > (tj  - t0)(Zj  - z0) 

Cd  »■(«'. P0»8,j  • (tj  - t#>* 

Observe  that  all  remarks  of  the  preceding  section  concerning  dT(pQ)  may 
also  apply  to  f(a',P0),  a situation  of  more  ramifications  than  are  at 
first  apparent.  In  addition,  it  might  be  observed  that  for  any  choice 
of  (Pj  | j * l,2,....n}  such  that  tj  - tQ  Is  constant  for  all  j among 
l,2,....n,  the  matrix  d y(a',p0)  becomes  of  rank  less  than  8,  and  In 
particular,  less  than  4,  hence  neither  of  the  Iterative  schemes  given 
by  (2)  or  (4)  will  converge.  Such  a situation  may  occur  In  any  of 
several  low  probability  situations,  depending  on  the  geometry  as  specified 
by  Uxj.yj.Zj)  | j * l,2,.,..n),  and  on  the  nature  of  (x0,y0,z0).  If,  for 
example,  (aj,a\a^)  * (0,0,0)  and  (Xj .y^ ,Zj)  are  points  on  the  clrcum- 
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ference  of  a circle  at  whose  focal  point  the  point  (x0.yo,zo^  Is 
located,  then  d .pQ)  Is  of  rank  3.  In  this  case,  however,  there 
Is  but  a single  point  of  discontinuity  near  (x0,y0,zQ).  It  Is  of 
Interest  to  note  that  In  this  case  the  condition  would  hold  wherein  the 
failure  of  a scheme  to  converge  could  specify  the  source  point. 

It  may  be  observed  that  If  an  extended  wave  generating  function  is  of 
the  type  specified  by  (5)  or  (8),  and  a suitably  restricted  estimate  of 
the  parameters  In  question  is  available,  the  solution  of  the  sound 
ranging  problem  Is  equivalent  to  the  transformation  d i'(a,,po)  being 
of  rank  8. 

In  practice,  restraints  on  the  physical  system  may  come  from  a primary 
estimation  scheme,  or  may  be  specified  by  the  user  in  the  form  of  a 

query  "Is  there,  based  on  knowledge  of  {p^l  j ~ 1,2 n)  and  a', 

a sound  source  near  p’2?"  If  so,  the  iterative  scheme  will  converge  to 
the  proper  value;  If  not,  the  scheme  will  diverge,  or  converge  to  the 
sound  source  Itself. 

CURRENT  INTERPRETIVE  TENDENCIES 

In  the  examination  of  behavior  characteristics  exhibited  by  the  proposed 
theory  In  a real  numerical  setting,  the  authors  decided  to  conform  with 
established  practice  and  consider  a six  sensor  system.  This  would  give 
rise  to  at  most  a system  of  six  nonlinear  equations.  The  reduced  system 
was  derived  from  a more  general  presentation  ("Sound  Ranging:  A Simple 
Model  and  Some  Ramifications")  by  Ignoring  vertical  variations  In  both 
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microphones,  sound  source,  and  wind  velocity.  The  variables  examined 
were  at  a minimum,  the  x and  y components  of  sound  source  position  and 
blast  time,  and,  at  a maximum,  these  variables,  together  with  wind  com- 
ponent velocities,  and  temperature. 

The  variety  of  possible  behavior  patterns  is  large,  even  in  such  an 
uncomplicated  sample  as  the  one  presented.  Differing  microphone  con- 
figurations, sound  source  locations,  meteorological  conditions,  and 
departures  in  preliminary  estimates  each  give  rise  to  differing  but 
overlapping  behavior  patterns.  More  than  a thousand  numerical  examples 
have  been  examined  for  the  purpose  of  clarifying  such  patterns,  and 
various  consistencies  have  begun  to  emerge.  Some  are  recognizable  at 
once  from  preliminary  analysis;  others  are  not  so  apparent.  These 
latter  results  are  useful  in  developing  a more  realistic  slant  toward 
further  analysis,  while  the  former  serve  to  verify  the  model’s  validity. 

In  the  particular  system  of  numerical  examples,  arrival  times  based  on 
sound  source  position,  microphone  positions,  and  wind  and  temperature 
information  were  calculated  to  the  limits  of  accuracy  of  the  Hewlett- 
Packard  9830A.  Among  the  Information  sought  were  possible  degree  of 
information  recovery,  number  of  Iterations  in  order  to  accomplish  re- 
covery, and  system  reaction  to  a progressive  degree  of  divergence 
between  "guess"  vectors  of  source  points  and  related  atmospheric  param- 
eters. Three  wind  velocity  situations  were  considered:  no  wind 

conditions,  winds  directed  toward  the  centroid  of  the  microphone  con- 
figuration, and  a limited  number  of  arbitrary  wlndspeeds  and  directions. 
Results  tend  to  agree  with  behavior  patterns  which  might  be  suspected. 
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but  some  interesting  phenomena  have  occurred  - the  action  of  which  is 
more  obscure,  and  which  will  be  of  considerable  interest  to  further 
exploratory  analysis.  A detailed  description  of  these  and  more  results 
will  be  presented  in  a more  visual  mode  in  a subsequent  report.  For 
the  present  situations  described  (i.e.,  hiqh  level  timing  accuracy), 
if  no  wind  was  encountered,  information  recovery  was  in  almost  all  cases 
complete,  even  when  extreme  departures  in  estimated  temperature,  wind 
velocity,  and  source  point  estimates  were  permitted.  In  general,  the 
number  of  iterations  to  recover  this  information  was  four.  Recovery 
was  considered  complete  when  estimated  and  known  quantities  differed  in 
the  hundredths  of  meters.  If  arbitrary  winds  were  introduced,  a mixing 
seemed  to  occur  between  temperature  behavior  and  wind  velocity  which 
indicated  that  error  behavior  and  stability  for  the  wind  and  no-wind 
situations  differ  considerably.  This  should  not  be  too  difficult  to 
verify  from  a purely  analytic  viewpoint,  and  may  possibly  require  an 
increase  in  accuracy  for  arrival  time  accuracies  when  wind  is  considered. 

The  greatest  problems  in  total  recovery  occurred  in  a purely  linear 
array  as  predicted.  The  various  other  arrays  which  were  examined  did 
not  appear  to  offer  any  really  significant  differences  in  behavior  in 
the  limited  cases  considered. 

The  documentation  of  the  above-mentioned  material  will  probably  be  fol- 
lowed by  at  least  two  reports:  a sensitivity  analysis  on  timing  and 

meteorological  errors,  and  an  examination  of  the  differing  behavior 
of  contemporary  sound  ranging  techniques  employing  data  derived  through 


simulation.  In  this  manner,  sufficient  statistical  control  is  available 
for  proper  data  interpretation. 

Note  that  each  sensor  configuration  requires  a certain  level  in  accuracy 
of  arrival  times  to  accomplish  recovery  either  of  a complete  parameter 
set  or  of  a source  point.  For  a given  configuration,  then,  knowledge  of 
system  behavior  may  be  specified  by  a set  of  "level  curves"  of  error 
values.  This  will  be  the  approach  taken  in  later  studies  when  the  effect 
of  configuration  on  error  behavior  is  examined.  In  addition,  more  gen- 
eral extended  wave  generating  functions  will  be  examined,  in  both  a 
determinative  and  random  setting.  When  all  ground  work  has  been  accom- 
plished, techniques  will  be  compared  by  use  of  a simulation  model  and 
then  with  actual  data.  This  should  place  the  techniques  as  described 
in  the  proper  functional  position  in  the  class  of  sound  ranging  schemes. 
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